Mapping the boundary of the first order finite temperature restoration of chiral 
symmetry in the (m w — m^)— plane with a linear sigma model 



T. Herpay* 

Department of Physics of Complex Systems, Eotvos University, H-1117 Budapest, Hungary 

A. Patk6s f 

Department of Atomic Physics, Eotvos University, H-1117 Budapest, Hungary 

Zs. Szep* 

Research Group for Statistical Physics of the Hungarian Academy of Sciences, H-1117 Budapest, Hungary 

o 

f^*) ■ P. Szepfalusy 5 

fN| \ Department of Physics of Complex Systems, Eotvos University, H-1117 Budapest, Hungary and 

Research Institute for Solid State Physics and Optics, 
Hungarian Academy of Sciences, H-1525 Budapest, Hungary 



3 



o 

(N 
> 



The phase diagram of the three-flavor QCD is mapped out in the low mass corner of the 
[m-K — mif)-plane with help of the SUl{3) x SUr(3) linear sigma model (LuM). A novel 
zero temperature parametrization is proposed for the mass dependence of the couplings away from 
the physical point based on the the three-flavor chiral perturbation theory (U(3) ChPT). One-loop 
thermodynamics is constructed by applying optimized perturbation theory. The unknown depen- 
dence of the scalar spectra on the pseudoscalar masses limitates the accuracy of the predictions. 
Results are compared to lattice data and similar investigations with other variants of effective chiral 
models. The critical value of the pion mass is below 65 MeV for all run values < 800 MeV. Along 
the diagonal m, = m K , we estimate m crit (diag) = 40 ± 20 MeV. 
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(JL, ■ 1. INTRODUCTION 

<D ' 

' The ambition of the exploration of the QCD phase structure corresponding to different breaking patterns of 
SUl{3) x SUr(3) chiral symmetry is the determination of the true ground state of the theory for an arbitrary 
set of quark masses m u ,md,m s in presence of a variety of intensive thermodynamical parameters, e.g. temperature 
(T), baryonic (hb), isospin (fij) and strangeness (fi s ) chemical potentials. The progress is continuous both in numer- 
■ ical lattice simulations [1] and in the application of effective models [2-5] for extracting results of phenomenological 
interest. The baryonic density of the Early Universe was very small when the cosmic expansion drove it through 
the stages of chiral symmetry breaking (the condensation of the different quark flavors) . Also for the extreme high 
energies of heavy ion collisions achieved at RHIC the average baryonic density of the final state is very close to zero. 
This motivates the present investigation where we concentrate on the case when all types of chemical potential vanish. 

Universality arguments [6] predict first order transition for m u = rrid = m s = and a second order one for 
m u = rrid = 0, m s = oo. One expects the existence of a triple point for some m s = m s , c . The most systematic effort 
seeking the explicit solution of the thermodynamics of the 3-flavor QCD is done with help of numerical simulations in 
the bulk of the (m u = m^, m s )-plane [7, 8]. However, by the nature of the lattice regularization, one explored to date 
mostly the region of rather massive u — d quarks, usually corresponding to pion masses of order 3-500 MeV (in these 
simulations m s is mostly kept fixed at its physical value). Lattice version of chiral perturbation theory (ChPT) is 
employed for extrapolating the results to the physical mass point. Also finite lattice spacing effects turned out rather 
important, therefore improved lattice actions gained significance in reaching physical conclusions. Common wisdom 
at present concludes that in the physical point temperature variations move the thermodynamical potential of the 
system analytically between the chirally symmetric and the broken symmetry regimes. 

At the same time constant interest is manifested concerning the location of the borderline of the region of first order 
transitions. If the border passes nearby, one might expect it to influence in a substantial way the transformation of 
the physical ground state [9, 10]. Numerical investigations were done and systematically improved for the 3-flavor 
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degenerate case m u = m d = m s ^ 0. The initial estimate of m crit (diag) « 290 MeV [11] was seen to be reduced to 
60 — 70 MeV [7] or may be to even further down [8] when finer lattices and improved lattice actions are used. 

Effective models (linear or nonlinear sigma models, Nambu-Jona-Lasinio model) represent another, in a sense 
complementary, approach to the study of the phase structure, which one expects to work the better the lighter quark 
masses are used [12, 13]. It is surprising that only moderate effort was invested to date to improve the pioneering 
studies of the SU(3) x £[7(3) linear sigma model by Meyer-Ortmanns and Schaefer [12] who used a saddle point 
approximation valid in the limit of infinite number of flavors, and derived m cr ; t (diag) < 51 MeV. An extension of 
their work to unequal pion and kaon masses was achieved by C. Schmidt [14]. He found m crit (diag) = 47 MeV and a 
phase boundary approaching the m^-axis rather sharply. The location of the tricritical point can be estimated from 
extrapolating his curve to rriK (tricrit) f» 70 MeV, although the expected power-scaling of the boundary curve with m r 
is difficult to disentangle from a simple linear regime. The phase boundary was calculated also by Lenaghan [15] using 
the Hartree-approximation to the effective potential derived in CJT-formalism [16]. For the complete determination 
of the couplings of the three-flavor chiral meson model he fixed the T = mass of the a particle in addition to the 
phenomenology of the pseudoscalar sector. The emerging phase boundary is rather sensitive to this mass. For instance 
in the case of 17a (1) anomaly, the tricritical kaon mass is n%K (tricrit) « 161 MeV (m s = 16 MeV) for m a = 800 MeV, 
and the expected scaling is not seen, while for m a = 900 MeV, mjf (tricrit) « 652 MeV (m s = 260 MeV). The 
estimate for m cr jt(diag) which one can extract from Fig. 3 of [15] for m a = 900 MeV is compatible with [12, 14]. 

In our opinion the greatest problem in refining the linear sigma model into a competitive tool of investigation of the 
chiral phase diagram is the difficulty of the determination of the quark (pseudoscalar meson) mass dependence of the 
couplings of the effective models. Almost all investigations tune exclusively the strength of explicit chiral symmetry 
breaking to cope with the variation of the pion and kaon masses via the Gell-Mann-Oakes-Renner relation. All other 
couplings are usually kept at the values determined in the physical point. One might note, however, some attempts 
to include also the variation of / w as deduced from lattice studies [13]. 

The novel feature of our paper is the parametrization of the couplings of the 3-flavor linear sigma model which 
ensures a full agreement with the results of ChPT for the variation of the tree-level pseudoscalar mass spectra as 
a function of the pion and kaon masses. In order to make the paper self-contained, we review in Section 2 the 
parametrization of the linear sigma model, which essentially follows Refs. [13, 17]. In Section 3 the relevant 0(1/ f 2 ) 
accurate results of ChPT [18-22] are summarized and used for the determination of the (m v ,mK)— dependence of 
the LerM-couplings. Full details of the parametrization can be reproduced with help of three Appendices. Next, we 
derive in Section 4 the equations of state for the nonstrange and strange condensates together with the gap equation 
for the common thermal mass which characterizes the finite temperature behavior of the scalar and pseudoscalar 
spectra. For this we use a variant of the Optimized Perturbation Theory [23]. In this way we partially avoid the 
imaginary mass problem of the standard loop expansions emphasized by [13]. In section 5 we argue that the phase 
boundary separating the region of first order transitions from the crossover regime varies sensitively depending on 
the assumption we make about the scalar sector when specifying the couplings of the model. Inspite of this variation 
we are able to conclude that the critical pion mass does not exceed 65 MeV in the region < tuk < 800 MeV. In 
particular the = mx diagonal is crossed by the phase boundary in the region 20 MeV < m cr ;t(diag) < 65 MeV. 



2. TREE LEVEL PARAMETRIZATION OF THE COUPLINGS 



The Lagrangian of the SUl(3) x SUr(3) symmetric linear sigma model with explicit symmetry breaking terms is 
given by [24] 

L(M) = ^Tr (dfiM^d^M + ^M f M) - /i (Tr (M f M)) 2 - / 2 Tr (AfUf ) 2 - g (det(M) + det(M f )) + e a + e 8 a 8 , (1) 

1 8 

where M is a complex 3x3 matrix, defined by the Oi scalar and m pseudoscalar fields M := —= ^(c; with 

v 2 i=0 

\i : i = 1 . . .8 the Gell-Mann matrices and Ao := y^fl- The last two terms of (1) break the symmetry explicitly, the 
possible isospin breaking term £3173 is not considered. 

A detailed analysis of the symmetry breaking patterns which might occur in the system described by this Lagrangian 
can be found in [13]. The fields <tq, erg both contain strange and nonstrange components. For the purpose of the 
exploration of the (m„ — mx )-dependence of the phase diagram we found more convenient to decompose the vacuum 
condensate into strange and nonstrange parts which is realized by an orthogonal transformation in the algebra basis 
and also defined the corresponding external fields: 
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ml = -^ 2 + 2(2f 1 +f 2 )x 2 +4f 1 y 2 + 2gy 


m\ a = -Mo + 2(2/i + 3/ 2 )x 2 + 4f L y 2 - 2gy 


m 2 K = -Mo + 2(2/! + f 2 )(x 2 + y 2 ) + 2f 2 y 2 - V2x(2f 2 y - g) 


ml -i£ + 2(2/! + f 2 )(x 2 + y 2 ) + 2f 2 y 2 + V2x(2 f 2 y - g) 


= ~d + 2(2/! + f 2 )x 2 + 4/ij/ 2 - 2gy 


= -Mo + 6(2/i + f 2 )x 2 + 4/iy 2 + 2gy 


m 2 vv =- f if ) +4f 1 x 2 +4(f 1 +f 2 )y 2 


ml m = Mo + 4/iz 2 + 12(/i + f 2 )y 2 


m L„ = ~ 2 9 x 


m L„ = 8/ix» + 2ga; 



TABLE I: The squared masses of the pseudoscalar nonet appear in the first column. The first two entries are the squared 
masses of pions and kaons, the last three rows represent the mixing in the r\ — rf sector. The second column contains the same 
quantities for the scalar parity partners. The phenomenological assignments of the scalar masses are discussed in the Particle 
Data Group (PDG) review on scalar mesons of Ref. [25]. 



where 



The fields with indices x, y appear in the matrix M as follows 



-V2 



M = 



75 X> 



V2 



V2 



diag(a x 



+ mt x , V2(a y + iiTy)). 



(3) 



(4) 



For the tree-level determination of the parameters of the system we have at our disposal the equations of state, the 
mass spectra of the pseudoscalar and scalar nonets and the consequences of Partially Conserved Axial- Vector Current 
(PCAC) relations for the weak decay of 7r and K . After some algebra (cf. [24]) one obtains the zeroth order term 



of the Lagrangian in the fluctuations around the expectation values < a x 
potential 



>= 



x , < a v >- 



y, which is the classical 



-L 



e x x - e y y - ^ + y 2 ) + gx'y + 2f lX Y + (A + 4fK + (/i + h)v 



Cy = V 



2 



The terms linear in the fluctuations must vanish, accordingly the two equations of state are 

dL 



E x := 



EL 



da x 
dL 



— X 

= y 



nix 



2gxy - 4f lX y 2 - 2(2 h + h)^ = 0, 
gx 2 - 4hx 2 y - 4(f 1 + f 2 )y 3 = . 



(5) 

(6) 
(7) 



The matrix of the squared masses can be read from the coefficients of the quadratic terms, see Table 2. There is a 
mixing in the x — y sector represented by entries of the last three rows. The mass matrix of rj fields is given also in the 
r/Q — -qs basis in Appendix A 1. The third and fourth order terms yield the three- and four-point interaction vertices. 
Finally, PCAC relates the condensates x and y to the pion (/„) and kaon (fx) decay constants: 



2V2f K = V2x + 2y, U 



(8) 



Equations (6), (7), (8) and those of Table 2 connect at tree-level the eight parameters of the Lagrangian (x, y, fj,o, fi, 
/2, g, e x , £y) and the physical characteristics of the meson sector, x and y belong to the coupling parameters of the 
shifted Lagrangian. The pseudoscalar masses and decay constants are better known than the corresponding quantities 
of the scalar sector, therefore the pseudoscalar sector is preferred over the scalars for fixing the parameters. The x, y 
condensates are simply obtained from (8). The couplings g and the combination M 2 := —(1q + 4fi(x 2 +y 2 ) can be 
determined by the knowledge of three pseudoscalar masses. The pion and kaon masses obviously should be selected 
since our purpose is to study the effect of their variation on the thermodynamics. For the third physical quantity, the 
trace of the mass matrix in the 77-sector is chosen, which will be denoted below by M 2 . 
This set of relations has the following explicit solution: 



X fir ) 

y = (2f K -U)/V2, 
(6/k - 3f n )m 2 K 



(2/x + U)ml - 2(f K - U)M 2 



Kf K -.U)(8f 2 K 

2f K m 2 K + 2(f K - U)ml 



8f K U + 3/2) 
(2/ K - U)M 2 



V2(8f 2 K ~8f K U + m 



(9) 
(10) 

(11) 
(12) 
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M 2 = ^M 2 + ^=(2f K - U)-2f 2 (U(U-2fK) + 2f 2 K )- (13) 

The above relations contain informations only from the pseudoscalar sector and were previously used in Ref. [13]. 

e x and e y are determined when using the above expressions in (6) and (7). But it is simpler to combine these 
equations in the Gell-Mann-Oakes-Renner (GMOR) relations and use the following equations instead of the equations 
of state: 

e x = mix, e y = — (m 2 K - m 2 n )x + m 2 K y. (14) 

These tree-level Ward-identities guarantee the Goldstone theorem at zero temperature. When to 2 = 0, the external 
field e x is zero and e y generates the nonzero value of rriK- The approach of taking into account the variation of the 
pion and kaon mass only by changing the external fields (cf. (14) ) was extensively followed in the recent literature, 
see e.g. [3, 12, 13]. 

The combination M 2 of /i and /x 2 is split up only in the expression of the admixed scalars, therefore the use of 
one characteristics of the mixed scalar spectra is unavoidable [17]. Nothing is known about the dependence of the a 
mass on m w and txik- We have applied the method described in detail in Section 4 also to the case when the mass 
of the a mass was fixed to a single value in the entire (m ff — itik )-plane. This scheme results in a phase diagram 
which is not compatible with the universal arguments on the nature of the phase transition in the chiral limit, at 
least for smaller sigma mass values, preferred nowadays [26]. Therefore some more flexible relation should be tested 
which allows the variation of the a mass with the pseudoscalar masses. We explored the consequences of assuming 
two different relations for the mass matrix of the scalars in the (m n — mif)-plane: 

Al. A first alternative is to assume that the mixing in the scalar x — y sector is absent (to 2 = 0), which along the 
m n = rriK line is the consequence of the U(3) x C/(3) Gell-Mann-Okubo (GMO) relation. 

A2. The SU(3) x SU(3) GMO mass formula for the scalars (A6) is fulfilled in the physical point with an accuracy 
of about —1.7%, supposing m a — 600 MeV. A second alternative is to require it to be fulfilled with the same 
accuracy for arbitrary m„,mK- 

Both assumptions involve a certain arbitrariness. The phase diagram was mapped out using both alternatives, and 
the resulting deviations give some feeling of the effects of our ignorance concerning the scalar sector. 
We give here the expression of /i and /i 2 , for the alternative ,A1' applied to the scalar sector: 

f{ Al) =-l y , , 2{A1) =Af^(x 2 +y 2 )-M 2 , (15) 

where the superscript ,Al' refers to the nonmixing of the x — y scalars. We can see in the equation above that in 
alternative ,A1' the coupling /i is directly proportional to the strength of the C/a(1) breaking determinant term in 
the Lagrangian (1). The implementation of the assumption ,A2' is more complicated, hence it is detailed in Appendix 
A2. 

The logics of the procedure sketched above can be summarized as follows: 

input: output: prediction: 




where m CT , toj are the mass eigenvalues of the admixed scalars and a is their mixing angle. 

The dependence of the parameters on and itik in (9) -(14) is not only explicit because one learns from the Chiral 
Perturbation Theory (ChPT) that all physical quantities (/„■, fx, M 2 ) featuring in this expressions also depend on 
m-ir ,rriK- Consequently, for the parametrization of the effective sigma model for arbitrary to t , ttik one should use 
the correct /u-(TO 7r , to^), fx (to^, ttik), M^m^^rriK) functions. In the next section we will construct these functions 
relying on results of the three-flavor ChPT. 
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3. DEPENDENCE OF THE COUPLINGS ON m w AND m K 



The fundamental problem of effective models in exploring the phase diagram of QCD in the (m ff — mx )-plane is 
the determination of the variation of the effective couplings when moving in the plane. The values determined in the 
physical point serve only as reference points, for a systematic exploration some reliable external reference is needed. 

The situation is somewhat analogous (but reciprocal!) to lattice QCD, where simulations are performed in a range 
of quark masses leading to much heavier pseudoscalars than in nature and some guidance is needed to arrive to the 
physical point. Chiral perturbation theory (ChPT) is used in this extrapolation [27]. Very recently it was applied 
in [28], for analyzing the pion mass dependence of the baryon masses of MILC collaboration. This suggests to us 
the idea to make use of ChPT results for deriving the parametrization of the linear sigma model away from the 
physical point. The issue of the compatibility of LaM and ChPT is not entirely settled. Recently the two models 
were compared in [29] in the light of the latest experimental data. The information available on the scalar sector, 
which improved considerably in the past few years, was used to fix some of the low energy constants of the ChPT 
with a more satisfactory result than thought possible previously. In this paper we fix the low energy constants within 
the ChPT, and adjust its renormalization scale, in order to match the pseudoscalar masses of the nonlinear sigma 
model with the tree-level spectra of LaM over an extended range of the (m OT — tuk )-plane. The calculated low energy 
constants of ChPT fall in the range commonly used in the literature. 

The essence of our approach can be understood by restricting our attention first to the functions / ff (m ff , itik) and 
fx {m w , ma) (the 77 — 77' mixing will be discussed afterwards). For this purpose, it is sufficient to choose the framework 
of SU(3) x SU(3) ChPT [18]. There 8 parameters (/, 4, q, Ma, L 4 , L5, Lq 7 Lg) were introduced, which determine 
m*,mx ) f.x,fK with 0(1/ f 2 ) accuracy: 



2A 



l + j2 1 



u 

Ik 



A(l+q) 



1 



H v + 164(2L 8 - L B ) + 164(2 + q)(2L 6 - L 4 ) 



+■ j5 \^Vh, + 84(1 + <z)(2L 8 - L 5 ) + 164(2 + q)(2L 6 - L 4 ) 



1 



1 + — (-2/^ -hk + 84L 5 + 84(2 + q)L 4 
1 ( 3 



1 + jz y- jGutt + fJ-n + 2hk) + 44(1 + q)L 5 + 84(2 + q)L 4 



(16) 
(17) 
(18) 
(19) 



where [ips — m 2 p s \n(mp s / Mq) / (32tt 2 ) are the so-called chiral logarithms at scale Mo, in which m 2 PS is substituted 
by the leading order expression for the squared mass of the corresponding member of the pseudoscalar octet. To this 
order one has in agreement with the Gell-Mann-Okubo formula m 2 = 24(1 + 2q)/3. It is worth to emphasize that 
Li do not vary with the pseudoscalar masses. 

The parameters 4 and q are related directly to the quark masses (cf. [18]) through 4 = B(m u + md)/2 and 
q = 2m s / (m u + mj) where B is determined by the condensate (uu) in the chiral limit. They can be expressed readily 
through the pseudoscalar masses and the chiral constants Li by 'inverting" Eqs. (16) and (17) to 0(1/ f 2 ) accuracy: 



4 = 



1 + q = 



l-jj\frr- l^r, + Sml(2L 8 - L 5 ) + 8(2m K + ml)(2L 6 - L 4 ) 



2m 2 K r 



1 



1 - -p (m>7 - + 8(77-4 _ rn 2 K )(2L s - L 5 )) 



(20) 
(21) 



It is sufficient to use the leading order relations of the two equations above to extract from Eqs. (18) and (19) the 
following 771^, TO/f-dependence for the pseudoscalar decay constants: 



fir = f 

Ik = f 



1 - -^(2/!^ + fijc - 4m^.(L 4 + L 5 ) - 8m 2 K L 4 ) 
f 

1 (Z 

1 - -p ( ^(Mtt + Mi) + 2fj, K ) - 4m^L 4 - Am 2 K (L 5 + 2L 4 



(22) 
(23) 



Using as input f n = 93 MeV, f K = 113 MeV, = 138 MeV, m K = 495.6 MeV, and m n = 547.8 MeV, and choosing 

M = AnU w 1168 MeV, / = 88 MeV, (24) 
one finds in the physical point the following values for the relevant chiral constants: 

L 4 = -0.7044 x 10" 3 , L 5 = 0.3708 x 10~ 3 , (25) 



which completes the continuation formulas for the decay constants (22) and (23). 
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These formulas enable us to predict the mass variation of the chiral condensates with help of Eqs. (9) and (10), 
and also the external fields e x ,e y from (14). The rax dependence of x, y, and e y is displayed for m^ = in Fig. 1. 
We remark that the only attempt, we are aware of, to take into account the nontrivial mass dependence of / w (m„.) in 
a thermal analysis, was based on fitting and extrapolating the mass dependence measured on lattice [13]. 



120 
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a? 60 : 
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FIG. 1: The tree-level kaon mass dependence of the T = condensates x and y and the external field e y for = {e x = 0). 

The chiral constants L§ and L$ are controlled by the values of A and q, respectively, taken in the physical point. 
Especially simple is the relation of L$ to the ratio q of the strange to average nonstrange quark mass. We take 
the value q — 24.8, which is close to the lattice determination and compatible with the range indicated by the 
PDG listing [25]: 20 < q < 34. For A we choose its leading order ChPT value in the physical point: A = A^-°\ Then 
using in the 0(1/ f 2 ) accurate expressions of A and q the phenomenological values of m 2 ,m 2 K with the Gell-Mann- 
Okubo formula for m 2 one obtains 



L 6 = -0.3915 x 10" 



L H = 0.511 x 10" 



(26) 



The values of the chiral constants Li, together with Mq and / can be used further for the continuation of A and q 
from the physical point to an arbitrary point of the (m ff — mx )-plane. 

The complete m v , m^-dependence of the couplings f2,g,M 2 given in Eqs. (11), (12), (13) requires also the knowl- 
edge of M^{m^,mK), that is the mass dependence in the (770, ?78)-sector, for which the application of U(3) x 17(3) 
ChPT is needed. The steps are quite analogous to what was described above, but the mass mixing makes it somewhat 
complicated. Since these formulae can be found dispersed in several papers we collect here the relevant formulae in 
more detail. 

In this sector, the 0(1/ f 2 ) ChPT results in a Lagrangian of the following form [19-22]: 



^os = ^zA.jd^d^rjj - - IK, 1 1, 11 , 



0,8, 



(27) 



where the elements of the real symmetric matrices A and D 



Ay — 6ij 



D..- D (°) 



(28) 



can be read off the papers [21, 22] and are compiled in Appendix B for the reader's convenience. The matrices a 
represent 0(1/ f 2 ) corrections to the zeroth order quantities. 

This Lagrangian is diagonalized in two steps. First one redefines the two-component vector rji as fj-, 
is followed by an appropriate rotation R(9 v )f]: 



A X J: 2 r\j which 



= R(0. 




= R(9 n ) • TO ^ 88 TO 'o8 



R- 



9 n ) = R(p n ) [I- -a 



D\l--a 



R- 



(29) 



Choosing 6 n 



ml , mz, , ml 

1100 ' I?08 ' V&S 



-20° and the experimental information on m 



one finds in the physical point the values of 



, which represent three relations restricting four chiral constants L 7 ,Vq 



(2) 



(2) (1) 



respective ChPT expressions for their masses. We choose the large N c relation v, 



(2) 



appearing in the 
-29. 3/ 2 [19] in order to have 
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as many unknown chiral constants as relations among them. This constant represents the contribution of the £/a(1) 
anomaly to the rj mass, dominantly determined by the topological features of the gluon configurations. It should be 
rather insensitive to the variation of the quark masses. From the expressions of the mass matrix elements listed in 
(B7)-(B9) one finds for the chiral constants: 

L 7 = -0.2272 x 1(T 3 , vP = 0.095, = -0.1382. (30) 

In the parametrization of LaM the sum of equations (B10) and (Bll) is used: 



1 



Ml = 2m 2 K -Zv {2) + 2(2m 2 K + ml){3v { 2 2) - V W) + -L 8v^ {2m z K + m*)(L 6 + 3L 4 ) + <G"„ - 3^) - imj,^ 



16 32 16 

+ y (6i 8 - 3i 5 + &L 7 )(ml - ra\ f + yI 6 K - 2< + n&m£ ) + —L 7 (ml + 2m\f 



(31) 



It can be checked that our results (22), (23) and (31) are the same as in [19], when the /Jps's and L4, L§, Lj, are 
set equal to zero (corresponding to the large N c limit). The chiral logarithm fi v contains the 77 mass at leading order: 

(m^° ) 2 = (Am 2 K — m 2 )/3, therefore the functions M 2 (m^^mK) and /^(m^, tuk), fidjn-K, ™>k) are only applicable 
when Am 2 K > rrQ. In addition we can rely on our "classical" approximation if the masses are lower than the chiral 
scale Mq. Eq. (31) together with (22) and (23) allows the computation of the couplings f2,g,M 2 in the pseudoscalar 
mass-plane. Their variation is illustrated in Fig. 2 for m, = 0. The theoretical quality of this parametrization is 
illustrated here by comparing m^m^m^ = 0) and 777^/(777,^,777^ = 0) as computed from the tree-level expressions 
of the linear sigma model with the results for the same quantities directly obtained from ChPT. Fig. 3 demonstrates 
that up to ttlk = 800 MeV the agreement is almost perfect. 
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FIG. 2: The tree-level kaon mass dependence of the param- 
eters of LcrM determined solely from the pseudoscalar sector: 
/2,g, and M 2 for = 0. 



PIG. 3: The tree-level kaon mass dependence of m v and m i 
for m-n- = 0. The labels refer to the results of ChPT and the 
predictions of linear sigma model (LcrM), respectively. 



The splitting of M 2 into f± and /Jq is realized in the mixing scalar sector, therefore it does not require any further 
consideration of ChPT. Their curves are shown for alternative ,Al' in Fig. 4, while the predicted masses of do and 
K are given in Fig. 5. It turns out that the alternative parametrization ,A2' leads to divergences in fx and /j§ for 
777a' — 777,^ ;$ 200 MeV. Therefore one cannot use it for the exploration of the whole (m^ — wlk) -plane. 

A final remark concerns the sensitivity of the T = mass spectra relative to the chiral constants (Li,v^ ). The 
values of the constants change considerably if, for instance, the large N c limiting formulas of ChPT are used. This 
change results in a rather large variation in the numerical values of fa, g, M 2 . However, the predicted masses of 77, r{ 
and the scalar sector remain almost unchanged. 

With this novel — mx -sensitive parametrization of the linear sigma model we are going to discuss the nature of 
the temperature driven chiral symmetry restoration in the following sections. 



4. QUASIPARTICLE THERMODYNAMICS OF THE SU(3) x 577(3) MODEL 



The aim of this section is to derive the equations of state (EoS) which determine the variation of the order parameters 
x and y with the temperature, including the existence of multiple solutions in certain temperature ranges. Results 
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FIG. 4: The tree-level kaon mass dependence of fi and (j.q 
with alternative ,A1' for the scalar sector, when m n = 0. 
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FIG. 5: The tree-level kaon mass dependence predicted for 
scalar meson masses ao and k in the linear sigma model, when 
= 0. 



of the numerical analysis of EoS determining the nature of the transition in function of the masses , m K will be 
mapped out in the next section. 

The renormalized EoS will be determined in the framework of Optimized Perturbation Theory of Chiku and Hatsuda 
[23] which starts by reshuffling the mass term of the Lagrangian density by introducing a temperature dependent 
effective mass parameter: 

L mass = ~M 2 (T)TrM^M + + M 2 (T))Tv MHl . (32) 

The first term on the right hand side is used in the thermal propagators of the different mesons. The second term 
in (32) represents the effective mass counterterm which is taken into account in higher orders of the perturbative 
calculations. 

The tree-level mass of it involves now the thermal mass parameter: 

ml = M 2 (T) + 2(2/! + f 2 )x 2 + 4f iy 2 + 2gy, (33) 

and all other meson masses to be used in the tadpole integrals below agree with the formulas appearing in Table 2 
with the replacement — yUg — > M 2 (T). If all quantum corrections are condensed into M 2 (T), then the tree-level masses 
of other mesons are expressible through the mass of the pion. One might expect that the pion has the lowest mass 
and therefore for M 2 {T) > these squared masses are all positive, which is not the case when — /Iq < figures in 
the propagators. We define a physical region of x and y where all tree-level mass squares are positive, and thus the 
one-loop contribution of the meson fluctuations to EoS is real. This region is most severely restricted by the masses of 
/o and a, which strongly decrease near the phase transition. We will look for the solution of the EoS's in the physical 
region. 

For the determination of the thermal mass we use the Schwinger-Dyson equation for the inverse pion propagator at 
zero external momentum. At one-loop it receives the contribution U(M(T),p = 0), which is the self-energy function 
of the pion at zero external momentum, plus the counterterm contribution — /x 2 , — M 2 {T). We apply the principle of 
minimal sensitivity (PMS) [23], that is we require that the pion mass be given by its tree-level expression: 

n w (M(T),p = 0) - (i 2 - M 2 (T) = 0. (34) 

H(M(T),p) itself is a linear combination of the tadpole and bubble diagrams (the latter not included in the treatment 
of [13]), with coefficients derived with help of the 4-point and 3-point couplings among mass eigenvalue fields. This 
step requires diagonalization in the (x,y) sector cf. Appendix C. The bubble contribution B(mi,m2,T,p = 0) at 
zero external momentum p can be expressed through tadpole integrals I (mi, T) as 

„/ m s /(mi , T) — /(ma, T) 

B(m 1 ,m 2 ,T,p = 0)= v u > 2 35 

m l — m 2 

therefore the self-energy is easiest to represent in form of a linear combination of tadpole integrals, which gives when 
substituted into Eq. (34): 

a— a. 7r 

= -M 2 (T)- f i 2 + J2 clJ(m m (T),T). (36) 

i— 7T, K, r},r} 
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Here c£ are the weights of the tadpole contributions evaluated with different mass eigenstate mesons a.; = o^, 7Tj . 
The integrals over the corresponding propagators are evaluated with effective tree-level masses where M 2 (T) replaces 
— /Uq. In this way (36) is actually a gap equation which determines the thermal mass parameter, M 2 (T). With help 
of Eq. (33) this equation can be also understood as a gap equation for the pion mass (the pion mass is present also 
in the expressions of I(m ai ,T) through m Qi !): 

a— cr, 7r 

ml = - $ + 2(2 f 1 + f 2 )x 2 +4f iy 2 + 2gy+ £ c^.I(m ai (T), T) . (37) 

i — 7T, K, Tj 7 T] 

Since this equation depends also on the order parameters x, y we have to solve in addition to the gap equation the 
two equations of state: 

-e x - fix + 2gxy + 4 W + 2(2 h + f 2 )x 3 + £™ , J l tlJ(m a% (T), T) = , (38) 



J i= 7T, /ST, 77,77 
K — (7,7T 

'i= 7T, it, 77,77 



^ y + ffx 2 +4/iX 2 y + 4(/i+/2)y 3 +E a ^ /(m ^ (T);T) =Q) (3g) 



with t*. and . giving the corresponding weights, listed in Appendix C. Ji is the isospin multiplicity factor: J v = 3, 
Jk = 4, and J VtV ' = 1. 

Equations (36), (38) and (39) represent a polynomial in x, y with divergent coefficients due to the divergences of the 
tadpole integral I. When compared to the expressions of the tree-level pion mass in Table 2 and the tree-level EoS (6), 
(7) one can uniquely absorb divergences into the couplings — /ig, /1, g. This step requires divergent counterterms 
as follows: 

. a (5/i + 3/2) A 2 (5f 1 +3f 2 )M 2 (T)-g\ A 2 

<V = 2 2 m "T2~' 

x Mh - h) , A 2 

69 = 2n 2 ln 7^' 

13/ 2 + 12A/ 2 + 3/ 2 A 2 

5/1 = 2^ ln 7^' 

c, 3A/2 + 3/I , A 2 
d/2 = 5 In To - , 

where A is the regularization cut-off and I is the renormalization scale. At T = with the replacement M 2 (T) — ► — /i 2 , 
these expressions agree with the known coupling renormalizations [24]. T-dependence appears only in the mass renor- 
malization, through M(T). Since they are proportional to higher powers of the couplings, this apparent environment 
dependence of the counterterm will be canceled by higher-loop contributions (see for instance, [30]). At the end 
of the renormalization we arrive at the same equations, just one has to replace /ig, /1, <?, I by their renormalized 
expressions (separate notation will be introduced below only for I —* I ft). 

The coefficients c£ . look at first sight horribly complicated since not only specific three-point couplings (see Appendix 
C) but also weighted factors proportional to (to 2 . — to 2 .)" 1 do contribute, cf. Eq. (35). However, a wonderful 
simplification occurs when working through this complicated expression, one finds cj. = Jit^./x. Then comparing 
the gap equation (37) to the EoS for the order parameter x, one recognizes the relation 

e x =ml(T)x(T), (40) 

which tells that the approximate solution constructed by us obeys Goldstone's theorem for the pions. This feature 
of the optimized perturbation theory was already emphasized in [23] in the context of the O(N) model. We mention, 
however, that when the symmetry breaking is realized by the appearance of two independent order parameters, the 
application of PMS in the form of Eq. (34) cannot keep the mass of the other pseudo-Goldstone boson, the kaon, at 
its tree-level expression. This means that the tree-level kaon mass does not satisfy the second relation of Eq. (14) and 
Goldstone's theorem. Had we chosen for the mass-resummation the self-consistent treatment of the kaon self-energy 
instead of pions, we would ensure that Goldstone's theorem is fulfilled for the kaons. Both relations in Eq. (14) can 
be fulfilled simultaneously only by resumming also one of the higher-point functions of the theory in addition to the 
mass. 

For the renormalization of I(rrii, T) we wish to use such a prescription, which allows to use further the parametriza- 
tion of the couplings realized with help of tree-level mass spectra. For this reason we decided to omit all temperature 
independent finite contributions from the tadpole and bubble integrals (the finite part of the 1-loop T = corrections 
to the self-energy), retaining in Ir only the contributions from the part of the integrands proportional to tib(uj,T), 
which is the Bose-Einstein distribution for a meson of energy to. The explicit form of the integral In(m, T) with this 
prescription is the following: 

1 f°° 

I R (m,T) = -^ dppn B (^/p 2 + m 2 /T). (41) 
^ Jo 
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Now one can proceed to the solution of Eqs. (37), (38), (39) for given m n ,mjc when T is varied. In the next section 
we describe in detail how first order phase transitions were detected and present the regions of the (m OT — )-plane 
where chiral symmetry restoring transitions take place with increasing temperature. 



5. THE PHASE DIAGRAM IN THE (m* - m K )-PLANE 



In this section we present our results on the phase diagram in the (m„ — m_R-)-plane paying a special attention to 
the physical point, the diagonal = tuk and the m v — axis. Investigating the nature of the phase transition 
along the diagonal is important because the result can be compared with lattice results [7, 11] and also with previous 
results [12-15], obtained in LaM. Moreover due to the degeneracy in the particle spectrum, the model is somewhat 
simpler on the diagonal, providing a good testing ground for our approximation. The m T = axis is relevant because 
of the the presence of the tricritical point which separates the region of first order phase transitions occurring for low 
values of mjf from the line of second order phase transitions. 

Since we have two order parameters: x (nonstrange) and y (strange), we have to monitor both of them in order 
to decide the nature of the phase transition. An interesting question arises whether one can speak about two phase 
transitions, one related to the melting of the nonstrange condensate and the other to the melting of the strange 
condensate. The sign for a first order transition is the appearance of three solutions for the equation of state (38) and 
(39) below a given temperature, corresponding to two minima and one maximum of the effective potential. 




FIG. 6: The temperature dependence in the physical point of: (a) the nonstrange (x) and strange (y) condensates; (b) the 
pseudoscalar (8 V ) and scalar (6 a ) mixing angles (in the (0-8) basis); (c) the mass of the chiral partners (ir, a) and (ao, 77); (d) the 
mass off fo,K,rj',K mesons. 



In Fig. 6 we present our results on the physical point, using alternative ,A2' with m a (T = 0) = 600 MeV. We 
preferred this one because alternative ,A1' gives m a (T = 0) f» 900 MeV, which is too high according to recent 
phenomenological studies [31] and experiments [26]. The evolution of both condensates at the physical point indicates 
a smooth crossover (see Fig. 6 (a)), with a peak in the susceptibility at around T = 210 MeV for the nonstrange and 
T = 310 MeV for the strange case. The evolution of the strange condensate is much slower. The restoration of the 
SU (2) x SU (2) symmetry can be seen by observing the degeneracy between the SU (2) chiral partners (tt, a) and 
(ao, if), Fig. 6 (c). We can observe the tendency of all the masses to converge at high temperature. Note, however, the 
gap between the two sets of chiral partners. This is the consequence of the J7a(1) breaking determinant term which 
enters with opposite sign in the expression of, for example, ir and ao masses. This is insignificant only for very small 
values of the strange condensate. The fact that, up to the temperature shown in the figure, the SU(3) chiral partners 
(it, ao) and (77, a) are not degenerate, indicates that the restoration of the chiral symmetry is not completed in the 
strange sector. We can also see in Fig. 6 (d) that the variation of the strange condensate is reflected the strongest in 
the mass of fo meson. 

The evolution of the condensates and masses is nicely reflected also by the temperature evolution of the mixing 
angles, Fig. 6 (b). The pseudoscalar and scalar mixing angles start at zero temperature at V = —10.45° and 9 a = 18.7° 
respectively, and they converge at high temperature. Up to the temperature we studied, they do not reach the ideal 
mixing ang le arcsin(l/ v / 3) ^ 35.264°, which means that fo and 77 are not purely strange mesons. In contrast to what 
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was obtained in [13, 32] the evolution of pseudoscalar mixing angle is nonmonotonic, it bends down and then up as 
the temperature increases. 

Next, we studied the phase boundary in the m n — m^-plane. As a reference, we considered the case when each of 
the zero temperature couplings of LuM has the fixed value calculated at the physical point irrespective of the value of 
rriir and ttik, except for the external fields e x , e y , which follow the variation of the and wijc according to Eq. (14). 
For m a = 900 MeV we obtained nearly the same phase boundary as in [15]. For m a = 600 MeV no phase boundary 
was found in [15] for > 0. With our method, the phase boundary is present, but it is not compatible with the 
universality requirement to have a first order transition in the neighborhood of the origin. In our view this represents 
an important argument for allowing the variation of all couplings with m T and m K . 
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FIG. 7: Phase boundary curves arising for LaM parametriza- 
tions compatible with T = ChPT. We present the boundary 
curves for alternatives ,Al' and ,A2' using the complete U(3) 
ChPT and the leading order large-iV c ChPT at 0{p 2 ). The 
exact phase boundary curve is expected to lie in the light-grey 
shaded region. Crossover (first order transition) takes places 
at the right (left ) of the shaded region. 
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FIG. 8: The temperature evolution of the nonstrange (x) and 
strange (y) condensates during a first order phase transition 
which takes place for m n — 10 MeV and tuk = 150 MeV. The 
curves were obtained for alternative ,A1' and the complete 
U(3) ChPT. Multivaluedness is observed in a given tempera- 
ture range in both condensates. 



Fig. 7 presents the phase diagram obtained in the case when all parameters are allowed to vary with and mx- 
Due to the uncertainties in the scalar sector and also due to different approaches of the chiral perturbation theory 
(three-flavor or large- N c ) we can give only a band indicative of the theoretical uncertainties concerning the location 
of the real phase boundary of the model. Note, however, that all variants give a first order transition near the chiral 
limit, that is for small values of both the pion and kaon masses. We see, that for any value of mx < 800MeV, 
the critical value of the pion mass does not exceed 65 MeV. Our estimate for the phase boundary on the diagonal 
is m cr ;t (diag) = 40 ± 20 MeV. In the figure also the line is displayed below which we cannot trust the results of the 
alternative ,A2', since its parameters diverge along the diagonal = mx- 

First order transitions are signalled by multivaluedness in the temperature evolution of both the nonstrange and 
strange condensates, see Fig. 8. For large values of the kaon mass, we claim that the phase transition is driven by 
the variation of the nonstrange condensate, since the apparently different solutions of the strange condensate are very 
close to each other, and all stay at high values. Subsequent decrease of the strange condensate at higher temperature 
displays only a crossover. Along the phase boundary the critical temperature reaches T c m 170 MeV near the physical 
kaon mass (in case ,A1' and using large N c ChPT), then drops to T c = 140 MeV at both ends of the mx range shown 
in the figure, which is due to the effect of the chiral logarithms. 

We could not provide evidence for a tricritical point on the m w = axis for any of the alternatives ,Al' and ,A2'. 
Alternative ,A1' seems to predict it for such a high value of the kaon mass, where one can not trust ChPT, while ,A2' 
does not work for = because the solution of EoSs leaves the physical region. 

Finally, we discuss a feature of our approximate solution in the low mass region which might be closely related to 
the problem of negative squared masses. It shows up the clearest along the diagonal, m w = mx, where the most 
plausible expectation would be to have a solution of EoS which satisfies, irrespective of the temperature, the condition 
o"s = 0. For the alternative ,A1', (,A2' does not work on the diagonal), it can be proved, using exclusively the tree- 
level stability criteria 3/i + /2 > 0, that going towards the origin below a certain value of the Goldstone mass there 
is always a temperature range in which one of the squared mass eigenstates in the mixing scalar sector has negative 
eigenvalue. In this range we find in the physical region only solutions with ag, ^ 0, that is the 'strange-nonstrange' 
symmetry is apparently broken in an intermediate temperature range. It is not clear if this phase corresponds to 
the absolute minimum of the free energy. This solution is characterized by a large difference between the mass mis- 
calculated from the second relation of Eq. (14) and the value of the pion mass, which is the largest one. In the mass 
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region where both solutions (with cr 8 = and og 7^ 0) exist the corresponding values of a are very close to each other. 
Therefore, we expect that even in the case where we cannot find the (true ?) minimum corresponding to erg = 0, a 
good estimate of the position of the phase boundary on the diagonal is provided by the as ^ solution. The problem 
of negative squared masses shows up also in the T — finite quantum correction coming from the tadpole integrals, 
which were omitted in this work. 

The above feature is a consequence of using tree-level expressions for the propagator masses. We certainly should 
have to go to higher-loop order in the resummed perturbation theory, also to take into account coupling resummations, 
for a complete resolution of the problem of negative mass squares including the assessment of the solution with og 7^ 0. 



6. CONCLUSIONS 



In this paper, we studied the phase boundary in the (m w — m^ )-plane allowing for the variation of all the parameters 
of the linear sigma model with and tuk- We used for this another low energy effective model, the chiral perturbation 
theory, which being a perturbative expansion in momenta and in quark masses about the chiral limit, provides, at 
each order of the momentum expansion, analytical relations displaying the dependence of the decay constants (f^, fx) 
and masses of the r\ and ?/ on the value of the pion and kaon masses. One could expect that the linear sigma model 
improved in this way will work reliably for small values of m n and itik ■ Using accurate formulas to continue from the 
physical point, this approach could become an alternative to the lattice which has difficulties in exploring this region 
when information would be available on the variation of the mass of the a or fa scalar mesons in the (m w — mx)- 
plane. The origin of the theoretical uncertainty of our findings is the lack of information on the scalar sector, which 
forces us to make assumptions. Lattice results about the mass dependence in the scalar sector would allow to reduce 
considerably the uncertainties of the parametrization of the model. 

The model was solved using a mass ressumation in the framework of the optimized perturbation theory in order to 
resolve the negative squared mass problem of the perturbation theory in the broken symmetry phase. Unfortunately, 
resumming only one parameter, the mass, while respecting the Goldstone's theorem for pions, violates Goldstone's 
theorem for kaons. It also does not solve fully satisfactorily the problem of negative mass squares in the whole 
mass-plane, since the absolute minimum might be located in the x — y-plane slightly outside the physical domain. 
Resummation of another coupling is needed to fulfill all requirements imposed by Goldstone's theorem. A possibility 
is to use the temperature variation of the coefficient of the Ua{^) violating term g. Motivated by lattice studies this 
possibility was investigated in [32]. 

Taking into account all theoretical uncertainties, we could estimate a band in the (to t — tuk )-plane for the phase 
boundary. Our estimate for the boundary point on the diagonal is m cr ;t(diag) = 40 ± 20 MeV, in nice agreement with 
the latest effective model and lattice studies. 



APPENDIX A: THE SU L (3) x SU R (3) LINEAR SIGMA MODEL AT TREE LEVEL 



1. Mass eigenvalues, and mass matrices in the 0-8 basis. 



Using the inverse of the transformation (3) , the mass matrix of ms can be written in the more conventional 770-778 
basis: 



<o = 5( 2 '<,+< H + 2 ^2"0 ( Al ) 

m ls = ( A3 ) 
The mass eigenvalues and mixing angle 8 r] are the following: 



"»W = IKoo + m lss T \/Koo - m ms) 2 + 4 <J' ( A4 ) 



2m 



2 

1708 



tan20„ = -5 22- (A5) 



where the ,-' sign refers to 77 and ,+' refers to r\ . These expressions hold also for the mixing in the scalar sector, 
where the lower mass eigenvalue is m^ and the higher is the squared mass of /o . 
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2. The ,A2' alternative 



For the scalar octet, there is a GMO mass relation similar to the pseudoscalar sector, which to leading order reads: 



4m K = m ao +3m a88 . 



We characterize its accuracy by the following quantity: 



4m 2 — m a 
3m L 



(A6) 



(A7) 



In the expression of r, the mass squares m 2 Q and m 2 are determined by /2, g, M 2 , which depend only on pseudoscalar 
mass squares: m 2 , m 2 K , and M 2 = to 2 + m 2 ,. In order to know nv 88 , we should have /i and /i 2 , separately. For 
this purpose, in the physical point, we choose an m^ h to get r ph . We require that the accuracy of scalar GMO to 
be independent of m-, m,K, that is r(m w ,rriK,'ma-) = rlrn^) =: r ph . After this, we can already determine m as& for 
arbitrary mx, m n from (A7), and we can split M 2 into /i and /Lt 2 ,. Using Table 2 and Eqs. (9)-(13) : 



fi 



(A2) 



rnMt^lk ± 104 ^/i ~ 58 /^ + 9 /> ~ 12 ( 2 /^ - ^Kfe - ^) 2 ) 
32(3/ 2 - 8Uf K + 8/|)(/a- - hf{r + 1) 

m 2 (r(16/, 3 + 8/f, + 24^/| - 39/ 2 / A -) + 4(/- + 2/ x )(/ J , - U) 2 ) 

32(3/ 2 - 8/-/* + Sf 2 K )(f K - / 7r )3(r + 1) 

M 2 {f K - U) 2 ((2f K - U)r + Ik - U) 
4(3/ 2 - 8.Uf K + 8f K )(f K - Uf(r + 1) ' 



# 2) (6/ 2 -8/-/ K + 8/ 2 )-M 2 , 
where r ph w -0.017, when m ph = 600 MeV. 



(A8) 
(A9) 



APPENDIX B: THE U(3) ChPT 



The elements of the mass matrix of 77-s are defined in (29) 

,2 



m 
m 



>/KS 

2 

>/(>o 
2 



^88^ + ^88 — (088^88 + a G9,D^) , 



Dqq + doo — (aoo-Do? + a 08-Do?) > 

^os + d 08 - kaosrf + _?&>) + D$(a 00 + a 88 )) , 



(Bl) 
(B2) 

(B3) 



together with the mixing angle defined in Eq. (A5): 



tan 29. 



2D 



(0) 



^00 -^88 



2<7, 



(18 



a m{D^ + -Dgg ; ) — -Dos ; ( a oo + a ss) , ^00 — ^88 — a oo-Doo + a 88-D 



)((!) 



(0) 
00 



,(0) 



2D, 



(0) 



D (o) _ D (o) 

^00 ^88 



(B4) 



The relevant expressions of the A, D matrices are given by [21], [22]: 

1 2 
pi' 
1 



aoo = -^-A(2 + q)(3L 4 + L 5 ), 



ass - -^(16^(2 + q)Li + —A(l + 2q)L 5 - 2ft K ), 



16 

y 



ao8 



16 1 

3V2T 2 



A(q-1)L 5 , 



D (o) _ 



00 

(0) 
88 

(0) 
08 



-3^ 2) + -A{2 + q), 



^A(l + 2q), 



-2\/_ 



A(q-1), 



(B5) 



^88 

doo 



1 2A 

_____ [QQA ((2 + q) 2 L 6 + 2(q - 1) 2 L 7 + (1 + 2g 2 )L 8 ) - (1 + 8q)n v - 9^ + 6/ijp] , 
-4A(2 + g)(4 1) -34 2) ) 

1 32A r 1 

A{(2 + g) 2 (6L 6 + 2L 7 ) + 6(2 + q 2 )L 8 ) - -fa + 6(1 + q)p K + (1 + 2^-) 



/ 2 9 



d 08 = A(q 1)2^4 1) - ^^A 2 (q 2 - 1)L 8 . 



(B6) 
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Therefore the elements of the mass matrix appearing in (B1)-(B3) are: 



1/88 



'too 



>/(>8 



32 32 

—A(l + 2q 2 )(2L s - i 5 ) - yi(l + 2?) (2 + q)L 4 



128 64 2 8 

-g-4(l - g) 2 L 7 + yA(2 + g ) 2 L 6 - -(1 + 8 g )/i„ - 2 M?r + -(1 + q)^ K 



(B7) 



-3vi 2) + ^A(2 + q) 



1 - 6?A 



18u 



(2 



1 



t~ (24^ 2) (3L 4 + i 5 



T 2 " 



y(2 + g ) 2 (3L 4 + i 5 -6i 6 -6i 7 ; 
yA(3(2 + q 2 )L 8 - (1 - <?) 2 £ 5 ) - -/Utt + 6(g + l)fi K + (2q + l)/i„) 



-2V2 



A(q-1) 



1 - 34 x) + ( 12t 



1 



(2) 



(B8) 
(B9) 



Carefully substituting the C(l// 2 ) accurate expressions of A and q from (20), (21) into (B7)-(B9) we find for the 
variation of m 2 oo , m 2 og , m 2 gg in the (m n — mA')-plane the following equations: 



4m|f — m 2 1 



'(88 



r- 



1 ^V., ., \™2 1 32. \/_2 2 \2 



32 

yieK - 2m* + mf^m 2 .) 
_ 3w ( 2 ) + 2< + < L _ 6v (D + 18u ( 2 ) + ^_ 24u (2) (i5 + ^ + 1_ 



-(-2/iK - /i^)m| 



- 7^)m 2 + y (2L 8 - L 5 )(m 2 - to 2 ,) 2 + + 2to^ 2 



16 

y 



2V2 
yr 
1 

7- 



(m 2 - m\) 



1-34 1) + -^(12^X 5 
1 



8(2L 8 -L 5 )(m 2 -m 2 ,)) 



4(L 4 - 4L 6 )(7« 7r + m K m n - 2m K ) + (-fi n -//*■ + fiK)m n + (-fi v - Mat) 



'A 



(BIO) 



(Bll) 



(B12) 



APPENDIX C: THE TADPOLE COEFFICIENTS IN EOS 

Below we list the nonzero three-point couplings, needed for the evaluation of the tadpole contributions to EoS (see 
Eqs. (38), (39)) in thex-y basis (4) 



TTyTTy 



'y u y 



2(2 h + f 2 )x 



1 a 

4/i2/ + 9 



K 2(2f 1 + f 2 )x-f 2 V2y+^g -f 2 y/2x + 4(/i + h)y 



2(2fr + f 2 )x 

4/ix 

-25 
2(2A + 3/ 2 ):r 



4/i2/ - g 
4(/i + h)y 


4/i y - 9 



k 2(2f x + f 2 )x + f 2 y/2y - j-g f 2 V2x + 4(/ x + / 



6(2/i +/ 2 )z 
4/ix 
8/i» + 2g 



H)V 



4/i y + .9 
12(/i + / a )y 
8/ix- 



(CI) 



The tadpole integrals are evaluated in the mass eigenbasis, therefore in the pseudoscalar x — y sector additional 
similarity transformations are needed in order to arrive at the coefficients of the r), tadpole integrals. The elements 
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of the 2x2 minor of t x ' y , which originally appears in the EoS, can be easily written in the mass eigenbasis. As an 
illustration we take the pseudoscalar (77, 77') sector: 



^XxG XX ~i~ tyyGyy -\~ t X yG X y Tr 



ta 



G X x G x y 



b yy , 



G x y G 



III) , 



= Tr 



R(6) i 



■ 2 
t 



R T (0) 



1111 , 



G(m v ,) 
v G(m n ) 



where G xx is the xx element of the 2x2 propagator matrix, R(9) is an orthogonal transformation defined by 
tan 29 = 2m^ xy / (m^ xx — ml ) which relates the x,y and 77', 77 basis. 

The diagonal (77,77), (77', 77 ) elements of the transformed matrix are the coefficients of the corresponding physical 
propagators and can be expressed as: 



1 



:(txx + tyy) i 



r]xx 



m riyy^(txx t yy ) + 2mf l xytxy 



(C2) 



4m 4 

rjxy 
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